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MTIONAL ADVISOEY COMMITTEE FOE AEEONAUTICS 
EESEAECH MSMOEANDUM 

VELOCITY DISTEIBUTIONS ON AEBITEAEY AIEFOILS 
IN CLOSED TUNNELS BY CONFORMAL MAPPING 
By H. E. Moses 


SUMMAEY 

Conformal mapping methods are applied to the calculation of 
the effect of channel (two-dimensional tunnel) walls on the ideal 
flow past arbitrary airfoils situated anywhere within the channel. 
The walls of the channel need not be plane but may have any shape. 
The results are compared in specific cases with those obtained by 
two approximate methods, of which the first is a first-order treat- 
ment using image vortices and doublets and the second is a higher- 
order correction developed by Goldstein. 


INTEODUCTION 

In reference 1 a conformal mapping method was developed whereby 
the ^zero-lift velocity distribution could be found for a symmetrical 
alrioil symmetrically located in a plane-walled channel. The purpose 
of the present paper is to extend the previous investigation to the 
case of an arbitrary airfoil situated anywhere within an arbitrarily 
shaped channel (two-dimensional tunnel). 

The Cartesian mapping function (CMF) , introduced in reference 2 
and used in the method of reference 1, is also used for the problem 
of the present paper. The velocity at any point on the airfoil in 
the channel is found in terms of the CMF and the known conformal 
transformation of a flat plate in a channel. The difference between 
this velocity and the velocity at the same point on the Isolated 
airfoil at the same angle of attack represents the effect of the 
channel walls. In order to obtain the velocity distribution on the 
airioil within the channel, the CMF is applied to doubly connected 
regions analogously to the manner in which Theodorsen's mapping 
function is applied in reference 3. 

The method is given, illustrated numerically by examples, and 
compared with corresponding results by the first-order image theory 
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and by the second-order image theory of Goldstein (reference 4) . 
In addition to the velocity on the airfoil, the veiocitj- on the 
channel walls is obtained by the conformal mapping method. 


SUffiOLS 


The more important symbols used in this paper are listed as 
follows : 


chord of airfoil 

section lift coefficient for isolated airfoil 
section lift coefficient for airfoil in channel 
distance between channel walls 
thicicness of airfoil 

undisturbed velocity at great distance from airfoil 

velocity on surface of airfoil in channel 

velocity on channel vralls 

velocity’' on isolated airfoil 

velocity correction, v^,^ - v^ 

angle of attack of airfoil 

angle of attack of flat plate 


h 

t 

T 

^cl 

Vc2 

^i 

Av. 

a 

a 

z 


effective angle of attack of airfoil with respect to 
curved streara 

plane of straight lines 

physical plane 

circle plane 
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METHOD OF COWFOBMAL MAPPING 
Th© CMP for One Contour 
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In previous applications of the conformal mapping method used 
in the present paper (for example, references 1 and 2), a single 
contour such as an airfoil in the physical plane (z-plane) was 
transformed into a single straight line contour (airfoil chord) in 
another plane (^ -plane). The line in the ^ -plane is related to a 
circle in a third plane, the p-plane, hy a knovm transformation 
that maps the unit circle with its center at the origin into the 
straight line such that the region outside the circle is mapped 
into the region outside the straight line. Because the contour in 
the z-plane also transforms into the same circle in the p-plane in 
such a manner that the regions exterior to the contours correspond, 
the function z - C _ is regular ©verjn^here on and outside the circle 
in tne p-plane. This vector difference z - ^ hetvreen conformally 
related points is called the Cai-tesian- raapping function (CMF)' 

The real and iraaginary parts of the CMF are denoted hy Ax 
and Ay, respectively. Because of the regularity of the CixlF outside 
the circle. 


z - 1 = Ax (p,qp) + iAy (p,<p) = ^ (l) 

where ^ 


P = pe 




^-n - "^-n + ^^-n 


On the circle p = e^*^ the folloiang relations hold; 

27t 

Ay (l,qp) = ^ I Ax ( 1 , 9 ') cot d9’ 

J 0 ^ 

-'2« f 

Ax (1,9) = - ^ I Ay (1,9') cot ^ d9' 

Equations (2) are the fundamental, equations whereby the transforma- 
tion between th© z- and ^-planes can be calculated. 
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The CMF for Two Contours 

In general, two contours in the physical plane, can he trans- 
formed into two straight-line contours in the t -plane. The lines . 
in the t -plane can, in turn, he transformed into two concentric 
circles in the p-plsme, whose centers are at the origin and whose 
radii are equal to 1 and q (q < l). The transformation is such 
that the region between the two circles is transformed into the 
region between the contours. 

In the case discussed in the present paper one of the contours 
in the z-plane is the airfoil itself; the other contour consists of 
the channel walls, both walls together being considered as one con- 
tour extending to infinity in two directions. The contours in the 
s-plane consist of a finite straight line into which the airfoil 
is transformed and a transformed channel whose walls are plane and 
parallel to the real axis. In the p-plane, the finite straight 
line, and hence the airfoil, are mapped into the outer circle whose 
radius is unity, and the channels of both the ^ - and the z-planes 
are mapped into the inner circle whose radius is q. Thus, as the 
outer circle is traced in a counterclockwise direction, the airfoil 
and the finite straight line are traced in a clockwise direction. 

In the same manner, as the inner circle is traced counterclockwise, 
the channel is traced clockwise. 

As in the case of the single contour, the regions at infinity 
in the z- and ^-planes correspond, but the vector difference z - t, 
is regular on the boundary of both circles and within the annulus 
formed by them. As before, z ~ t is the CMF and Ax and Ay are 
its real and Imaginary parts. Because of the regularity of z - t 
in the annulus, the CMF may be expanded as follows (cf. equation (1)) 

00 

z - t = Ax (p,cp) + iAy (p,cp) = C^p^ (3) 

where 



Inasmuch as the full Laurent series is used in equation (3), the 
relations between Ax and Ay on the two circles differ from the 
simple relations given by equations (2). 

Appendix C of reference 2 provides relations between the com- 
ponents of the CMF on the two circles, but the expressions are not 
easily used for the purpose of calculation. More convenient rela- 
tions have been derived in reference 3. Although the correct result 
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is obtained, the method of derivation is not fully given. The rela- 
tions aci’e derived in more detail in appendix A of the present paper. 
These relations between the components of the CMF are the following. 


The subscripts 1 and 2 indicate the values of the CMF on the 
circle of unit radius and the circle of radius q, respectively. 
That is : 

AXi(cp) = Ax(l,cp) 

Ayi(<¥>) = Ay(l,cp) 

AXgCcp) = Ax(<i,cp) 

AygCcp) = Ay(q.,cp) 

Then, as shown in appendix A; 

•2xt 
0 


(4) 


,2« 

Ax3_(cp)= aQ+- J Ay3_(cp') 


1 cot , 

2 2 ^ 


\'2a^“ 

> sin n(cp' -cp) 

1-q 


dcp'l 




dcp' 


1 

AxgCcp) = ao-- J Ay„(cp') 


1 - -u (<P'-'P) \ 2q2n . , . ^ 

_ cot siti n(cp'-<p) 


dcp' 


^ ,2jt ^ n 

- I Ay,(cp')\ — 3^ sin n(cp' -cp) dcp' 

” ■'o 1-q^" 


‘(5) 


1 f 

Ayi(cp) =1^0 - - j AXj^(cp' ) 


0 

2it 


1 ot ?afl 


.n 


if ax2(cp.) 

0 l-q 


2 ' / T 2n 

1-q 


sin n(cp' -cp) dcp' 


sin n (cp ' -cp) 


dcp' 


1 / 

Ay2(cp) =bQ + - J AXgCcp') 


0 

2n 


1 --4- ('P'-cp) \ 

-OOt 


sin n(cp ' -cp) 


dcp' 


- - / Ax, (cp' sin n(cp'-cp) dcp' 

” 0 Z_.l-q2n 
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and 


rzti j2n I 

Axi(cp) = AX2(cp) dcp i 

Jo Jo ! 

> ( 6 ) 

f2rt f2jt 

<3;^ = j A72 (<■■?) <3cp | 

./o Jo 


The introduction of elliptic functions simplifies equations (S). 
The elliptic functions introduced at this point and used at othoi’ . 
places in this, paper are treated in various texts with varying 
notations. The notation used througiaout this paper is that of 
Tannery and Mollc (reference 5). From reference 5 (t. IV, p. 100 ), 
the following series for the i5-f unctions are obtained: 


i = i cot TtV +^y ^ ■ sin 2itnu 

2rt 2 ^ l-q2n 

^ - 1 


i_ JfL 

2rt 2n 

1 


sin 2nnu 


(7) 


( 8 ) 
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Hence, from equation (s) : 


Axi(qp) = ao + -i- 


rzn 


0 


Ay-i ( 9') 


/ :p' 

,x I j 


V 2n y 


dp' 


r2rt 




— P j Cpo — - -— ■■■ ' ■ d 'P' 

2jt2 Jo ^ 




(£L2\ 

V 2n / 


AxgC'P) = Oq 


f2rt 


2rt2 jo 


AygCcp') 




\ Zn 


l 2n / 


- ^ V’ 


f 2n 


2 l( 


2 n^ Jo 




h' 


(J£LJ£) 

\ 2n / 




!£sl\ 

'4 \ 2jt / 


d 9 


Ayi(9) = to - 


2«2 J 


{^) 

Ax^Ccp') -1 — l_2£_l 


d'C?' 


(9) 


r2jr ^ 

I AXgC-p’) -- - ' d(p' 

2«2jo 2 

^ V 2 rt i 


f 2it 


AjoC'?) = to + I AX2((p’) 
2^2 Jo ^ 


a ■ (^] 

,„.) t vaW arp' 


’’i 


2jt" 


'2n 


, / cp'-9 


Ax-, ( cp') 


-^4’ I 2:t 
4 V 2rt y 


d <P' 
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In the same \ra,y that the relations expressed hy eq.uations (2) 
are a limiting form of Poisson's integral, the integrals in 
equations (5) or (9) are limiting forms of Villat's analog to 
Poisson's integral. Villat's integral (reference 6) gives the value 
of a function within an annulus when the real part of the function 
is Icnown on the hounding circles. 

The relations expressed by equations (5) reduce to those 
expressed hy equations (2) when the radius of the inner circle 
approaches zero; that is, when the channel walls move to infinity. 
The signs differ, however, because the CiiF that is defined within 
the annulus in equations (S) is defined within the outer circle as 
the radius of the inner circle goes to zero, whereas, in the case 
of equations (2), the CMF is defined outside that circle. 


The ^ -Plane and Its Transformation into the p-Plane 

As already described, the ^ -plane contains a plane-walled 
channel within which there is a flat plats. The transformation ■ 
mapping these contours into two concentric circles has teen obtained 
by Tomotika (reference 7), who has also obtained the velocity poten- 
tial for this case, Tomotiica's results will be briefly presented 
and the form in which they are most useful in applying the CMF method 
will be given in more detail. 

Let a be the angle of attack of the flat plate. The trans- 
formation between the ^-plane and the p-plane is shown in figure 1 
and given mathematically as 

Cc 

f = - £il ,e \2 / \ 

Z_n(l + 2q2^ cos 25 + q^=") 

1 

X [p^(e^^ + q^^e“^^) + p"^(e"^^ + q^^e^®)| + T (lO) 


or in another form 



( 11 ) 
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where h is the distance between the channel walls, 6 = ^ - a, 
and T is a constant. The substitution of p = e^'^ in equa- 
tion (10) yields the equation of the flat plate; the substitution 
of p = qe^^^* in equation (ll) yields the equation of the channel 
walls, which are parallel to the real axis. The use of the two 
forms of the transfonaation simplifies the resulting equations 
in ip for the flat plate and channel walls. 


Four values of the central angle p (p^, P 2 > andp^) are 
important in the mapping. From, reference 7 the points p = 

i (2jt-pi ) 

e , denoted by B and B* map into the stagnation points 


on the flat plate for zero circulation; the points 

i (Sn-pp) 

qe , denoted by H and H' map into tree , 

ip^ IPa 

points p = e e denoted by A' and A, 

extremities of the plate. The points are shown in 



respectively; 

map into the 
figure 1. 


the 


The values of q and p 2 in equations (9) and (10) are 
determined by the length of the plate, its position, and the various 
relations between the four special values of p . From equation (lO) 
or equation (ll) " 


P 3 +P 4 = 2 a 

( 12 ) 

/Pl-*P 2 \ /Pl-P 2 \ 

»4’V2n ) 5 

(13) 

4\ 2 it / 2n J 


00 

V 

^ y' ^ 2 n ® 



1 


h\-^) = %\-27-y 


(3.4) 


L 

h 



n(p3~p4) 

q^ sinnPg sin g 

n(l-2q^^ cos 2a + q^^) 


fcos(n-l) a -q^^ 


cos(n+l) a] 


1 


( 15 ) 
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^ = i / 21 - q>^\ + — sin a / 
n n \ 2 It / 


„ . n(93‘'?4) 

q.^ sin ncf^ cos g 


la + 


/n(i- 2 q^^^ cos 2 a + q‘ 


X 


j sin(n-l) a - sin(n+l)^ 


( 16 ) 


vihere L is the length of the plate and d is the vertical distance 
of the midpoint from the center line of the channel. 

In order to find the values of q, 
a given length and position of the plate, equations ( 12 ) to (16) 
should he solved simultaneously. 


In principle, it is possihlo to transform the. flat plate at 
any value of a into the airfoil at angle of attack set at a. 

The value of a is fixed at the value that accomplishes the trans- 
formation with the least labor. In the case to he calculated, 
a is set equal to zero. For this value of a, Tomotika's formulas 
(reference 7) are considerably simplified. If the distance between 
the channel walls is taken as unity, the equations simplify as 
follows : 

(pn.p-n, , , . (17) 

“ na-l 2 “) 

or equivalently 


Cx> 


r _ ^ !_±. 

s “ -rf \ o_. 


icpo 

1 , _ P-ae ^ 


\ 

,/ 


sin nf% |7pf /pV”|, , , 

w ’V5y Jf 


/ n(l-q^^) 

1 




(18) 


and Cpg- found by using the equation 


d 



(19) 


The quantities q, 93 ^ and T 4 may he found hy solving 
simultaneously equations (12)^ (13); (K); and (15) , which also 
become simpler than the equations for the general case. The 
constant t may he restricted to real values; because it merely 
determines the position of the channel and the flat plate ^nth 
respect to the axes in the ^-plane. 
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A special case useful for numerical work is that for which 
d = 0. For this position of the flat plate = Pg = P3 = -^4 = Jt/2 


L 


„2n-l 

"A. (2n-l)[l-q2(3n-l)j 


( 20 ) 


Henc0^ q can be found froin ecjuation (20) alone when the lengt^^^ 
of the plate is prescribed. 


^ If ^ is separated into its real and imaginary narts^ 

C = ^ + iT)^ ^ and t) can be found as functions of 9 . The 

e(iuation of the flat plate is found by setting p = e^'^ in equa- 
tion (17). Then^ when the subscripts 1 and 2 denote the' values of 
the function on the plate and on the channel^ respectively^ 

00 

q^sin nqjg 

1 

^1 = ° ( 22 ) 


Thus the flat plate lies on the real axis of the ^-ulane. The 

equation of the channel walls is found by setting p = qe^'-P in 
equation (18). Then 


io = 


? i 


sin I (cp-cpg) 


sin (rp+cpg) 


Z q^^sin nQ2 . 

4 > Sin n'P 

n(l-q2n) 


1 


+ T (23) 



^2 



1 = 


/I 

V2 



(2n-P2 >^P > '{>z) 


> 

I 


( 24 ) 
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Airfoil Position and Adjustiaents in Terms of the Q'®’ 

The z-plane and the ^ -plane are shown superimposed in figure 2 
in which the geometric meaning of the CMF is also indicated. If 
the abscissas and ordinates of the airfoil are denoted by 
and the abscissas and ordinates of the channel walls by yg; 

the definition of the CMI’ shows that 


2ii(cp) = ^2^(9) + AX3 _(<p) 
y3_(cp) = Ayj_(cp) 

X2(cp) = igC?) + AXgCc?) 
yg(9) =^ 12 + Avg(9) 


(25) 


In order uo determine the constants q. and ’"P 2 that appear 
explicitly in the expressions for t and tj and also the angles Cpjj 
and that correspond to the leading and trailing edges of the 

airfoil, the airfoil is placed in a normal position with respect to 
the y axis. If c is the chord of the airfoil and a ;s the angle 
of attack, the normal position is given by 


- — cos a 
2 


— cos a 


(26) 


From equations (25) and (26), the following formula is obtained: 

t ]_('■%) - ^i('Pw) = c cos a - Axp(9t) + Axp(cpff) (27) 


The angles cpjj and 9,j corresponding to leading and trailing edges 
are obtained from the condition of a maximum for the abscissa x-| (cp ) 


dXj (ipi'j) ^ ^ 

dcp 

or, by equations (20) and (24), 


dxj^(cp^) 

d9 


= 0 


( 28 ) 


dAX]_(cp) 

dcp 


2 

4 \ q’^ sin ncpg 
1 - q2n 

1 


COS ncp 


(29) 


for 9 = 9 jj or 9tp. 
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The value of cpo, or what is equivalent the value of d, is 
found from 



C 2jt 


Ayi(<p) = i 
Jo 


Ayg(^) ^<9 


( 6 ) 


as follows: Let r(<p) denote the value of the ordinate of the air- 

foil measured from the center line of the channel in the ^ -plane. 
From the definition 


r(cp) = Ayj_(cp) + d 


(30) 


Hence^ using equation (6)^ 

, f 2jt , r 2:t 

A = ^ i-6p) dC? - ^ I dcp (31) 

Jo Jo ^ 

and cp2 is obtained from equation (IS). 

The constant t is obtained by adding the equations of (26) . 
The resulting formula is 

CO 

4 \ sin n(fe n(cpjT+cpp) n(cpjr-cpg,) Ax, (cpj + Axn(q)rp) 

T = - > sin — i_ cos i ^ 

^ , n(l-q^^) 2 2 2 

1 (32)' 

These equations completely determine the constants q, T, 
cpU, '92 terms of the CMF. The value of 'Pg calcu- 

lated from equations (3l) and (19); the use of this value in 
equations (27) and (29) permit these equations to be solved simul- 
taneously for q, cpjj, cp^; and finally t can be determined from 
equation (32). 


Velocity Distribution on the Airfoil 
and on the Channel Walls 

The complex velocity potential W^ derived from the results 
of reference 7, is 
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where V is the velocity at infinity and P is the circulation. 

The velocity distrihution on the airfoil and on the channel 
walls is obtained from the velocity potential given by eq.uation (33) . 
The formula for the velocity in the z-plane is 


= 


1^ 

dp 

dz 

dp 


On the airfoil, from equations (8) and (33) 


(34) 


f§M\ 


le 


-i<p 


d'-P 


2rt 


/ sinncpp^ \ 

j QY y — — ^ cos ncp> - r 




7 


/ 


(35) 


The circulation P in equation (35) is adjusted to satisfy the 
Kutta condition at the trailing edge of the airfoil = 0. 

The result is •■P=iPrji 


00 


q^ sinn rpg 


P = 8V 


1-q 


2n 


- cos n^piji 


(36) 


Also 


/ii\ = 

(dp) Vdp/ 


■D=e ^ 




- 1 


^-i,p dAy;^(9) 


dip 


+ e 

dcp 


Oo 


n 


41e-W \ 4“ _ IdA^ 




1-q 


2n 


d-p 


dcp 
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Hence, the velocity distribution on the airfoil is 

00 

sin ncpg 


cl 4 


1-q: 

1 


2n 


(cos ncp - cosn'^ip) 


(38) 


-2 


/ |4 'X <1^ sin ncpg 


/ , 2n 

/ i / — •• 1-qi 

/ L 1 

where v^,-) has been witten for v^. 


dAX]_(cp)j I dAy^('qp) 




drp 


+ ! 




L 


The velocity distribution on the channel walls is found by 

replacing p by qe^*^ in equation (33) . The substitution results 
in 


m 

\ap/ 


/ 

* ■ 


iq-^e 


sin Cp2 


,2n 


p=qe 


iCp 


cos 9-COS 1^2 
\ ^ 


4 x^ sin ncpg cos nCp p 


, \ 


l-q' 


2n 


2nYj 


( 39 ) 


where P has the value given by equation (36). 



dAxgC-p) _T dAygC'?) 

. . -u n — 
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Lift on the Airfoil in the Channel 

The lift on the airfoil in tiie channel can he found hy evalu- 
ating a Eiodified form of Blaaius' integral in the p-plane. The 
expression for tho lift involves the CMF and the radius of the inner 
circle; that is, the lift depends on the shape and the position of 
the airfoil and on the shape of the channel walls as well as on the 
circulation. This dependence is in contrast to the case of the 
isolated airfoil, in which the lift on any body is the same for a 
fixed circulation. The dependence of lift upon the airfoil shape 
for the case of the airfoil in a plane -walled charmel has also been 
shoT-m by Havelock (reference 8) who finds the potential function 
directly without the use of conformal mapping. 

The expression for the lift is too complicated for numerical 
calculation. A more convenient way of obtaining the lift is to 
integrate the pressure distribution on the airfoil or the pressure 
distribution on the walls. 


Method of Successive Approximations for Obtaining CMP 

The CMP can now be calculated for a given configuration by a 
method of successive approximation analogous to that of reference 2. 

1. The airfoil and the channel walls are drawn such that the 
airfoil is in the normal position, as shown in figure 2. The center 
line of the channel in the ^-plane is located on the figure in order 
that the airfoil 'ordinates r(<p) may be read. The scale is so 
chosen that the distance between the channel walls in the ^-plane 

is unity. 

2. Prom a previous approximation, approximate values of q, 

T, cp 2 ; and are imown, as well as approximate values of the 

abscissas x-j_(cp) and X 2 (CP) at a convenient set of values of cp 
from 0 to 2n radians. Through the use of the known values of XQ_(cp), 
r(<p) is measured. A set of values of Ay.,(‘P) are also measured 

through the use of the Icnown values of X 2 (<P). A value of d and 
a new value of cpg are obtained from equations (31) and (19) . 

If no better values are available, the initial approximation 
for Xj_(cp) and '>^{'■ 9 ) may be that obtained for the flat plate 

situated along the center line of a plane-walled channel. In this 
case Xj^ and X 2 are given by equations (21) and (23) for 
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and ^2* value of q. is obtained from equation (20), where L 

is replaced by c cos o.» Both CP 2 and cpjj equal Jt/2 and 
equals 3n/2. The constant T equals zero. 

3. The functions and AX 2 are calculated by means of 

the first and second equations of (5) . The value of q used is 
the approximate value of step 2. The numerical details of the 
calcu?Lation are given in appendix B. 

4. New values of cPj.j, and q are obtained by solving 

equations (27) and (29) simultaneously for these quantities. 

An alternative method of determining and q is a 

purely graphical one. The approximate function X;j^('p), which is 
also a function of q, is plotted against <p in the regions of 
the extreme values of x-j^. Fi’om this graph 'Pjj and spip are 

detemined. These values are substituted in equation (27), from 
which a new value of q is obtained that is used to re-evaluate Xj^. 
The procedure is continued until sufficient accuracy is obtained. 
Finally T is calculated from equation (32). 

5. A new set of values for xi(<p) and X 2 (q)) a^re calculated 
using the new values of the constants and the values of Ax^^ and AX 2 
calculated in step 3. 

Steps 2 through 5 are repeated until a plot of y('?) against 
x(cp) for both the airfoil and the channel walls yield shapes that 
are as close as desired to the sliapes plotted in step 1. 

If the walls of the channel in the z-plane are flat, Ay 2 (<P) 
is sot equal to zero, and a considerable simplification in the 
numerical procedure results. This case is the most common and the 
method is not at all difficult to apply. The discussion of numerical 
results will provide an idea of the actual work involved. 

After the components of the Ct-IF and the various constants have 
been evaluated by the method of iteration just described, the veloc- 
ity distribution may be found from equations (38) and (4l) for the 
airfoil and for the channel walls, respectively. The derivatives 
of the CMF in. the formulas for the velocity distribution were 
measured in the cases calculated j although an expression exists 
that gives the values of the derivative in terms of the CMF as in 
reference 1, it is too cumbersome to use. 
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ILLUSTRATIVE EXAMPLES USING C0NF0Rf4AL MAPPING 

The method, of conformal mapping outlined has been applied to 
the 12 -percent symmetrical airfoil treated in reference 1. The 
ordinates of this aii-foil are given in table 1 and the airfoil shape 
is shovn in the figures in which the velocity distributions are 
plotted. For the calculations of the present paper the. airfoil was 
assumed to be placed at the center of a plane-walled channel. The 
chord to height (c/h) ratio was taken to be 0.5. Velocity correc- 
tions wei-e calculated for angles of attack of 0° and 4°. 

For the case of a = 0° the range of from 0 to 2h radians 
was divided into 24 equal intervals. T\to approximations, stai'-ting 
from the x(<p) of the flat plate, wore necessary for the derived 
airfoil contour to coincide with the given contour for a scale of 
chord length of 20 inches and ordinate scale five times that of the 
abscissa scale. In no case were more than six tei-ms used in any of 
the infinite series in the preceding formulas, for the series con- 
verge rapidly. The velocity distribution for the case of a = 0'^ 
is shown in figure 3. The velocity distribution on the walls of the 
channel is included in the figiure and is drara to a scale five times 
as large as the scale for the velocity distribution on the airfoil. 
The Cf'lF together with the velocity distribution is given in table 2. 
The velocity distribution on the airfoil for this case had been 
previously calculated by the method of finite chord in reference 1, 
The results are compared in figure 4 and are in close agreement, 
which indicates that the numerical methods used in both processes 
were sufficiently accurate. 

The velocity distribution for the case of angle of attack of 4° 
is plotted in figure 5. Figure 6 shows for the purpose of compai-ison 
the velccit;/ distribution for the airfoil in the free stream at 
a = 4“^. In this case four approximations, starting from the flat 
'plate, wore necessary to obtain coincidence between the derived air- 
foil and the given airfoil to the same ordinate and abscissa scale 
as in the case of a = 0°. In the first three approximations the 
9 range was divided into 24 equal intervals, but in the fourth 
approximation the length of the intervals was halved so that the 
CMF was evaluated at 48 points. The mapping data and velocity dis- 
tribution are given in table 3; the, nature of the C1>1F is shov/n by 
figure 7 where the component fiuictions are plotted. The velocity 
distribution for the airfoil in the free stream was obtained by the 
method of reference 2. ' 

The velocity correction for the airfoil at an angle of attack 
of 0° was discussed in reference 1. The velocity corrections for 
the airfoil at the angle of attack of 4° are plotted in figure 8. 
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The irregularities of the correction are due to local curvature 
fluctuations of the airfoil surface and correspond to the irregu- 
larities found in the corrections for the same airfoil at a = 0°. 

( See reference 1 . ) 

The velocity corrections are positive on the upper surface of 
the airfoil hut are for the most part very nearly zero on the lover 
surface. This behavior of the correction indicates that the lift 
on the airfoil in the channel is greater than that on the airfoil 
in the free stream.. The increase in lift has been shovn by other 
authors thi’ough the use of approximte methods (see references 4, 

1 , 8f and 9) and will be further discussed. 

The influence of the airfoil on the velocity distribution on 
the channel walls is shown in figures 3 and 5. The velocity dis- 
tribution on the walls is very sensitive to the angle of attack. 

\^/hen the angle of attack is 0° (fig. 3) the nondimens ional velocity 
on both the walls is greater than unity. The velocity rapidly 
approaches unity both upstream and downstream of the airfoil until 
at 1.75 chord lengths upstream and downstream of the origin the 
velocity has decreased from its maximum value 1.03 to substantially 
the value 1.0. 

In contrast, when the angle of attack is 4° (fig. 5), the 
velocity is less than unity on the lower wall, and on the upper wall 
the velocity markedly increases over the velocity for the case of 
a = 0*^. The mximum velocity on the upper wall moves forward toward 
the position at which the airfoil approaches closest to the wall; at 
the same time the minimum value on the lower wall is located at the 
position near the leading edge where the zero streamline rises to 
meet the airfoil at the stagnation point. On both the upper and 
lower walls the velocity approaches unity less rapidly than in the 
case of a = 0*^. On the upper wall the maximum velocity is 1,035; the 
velocity 1.75 chord lengths upstream of the origin is 1.013; the 
velocity 1.75 chord lengths downstream is 1.010. On the lover wall 
the minimum velocity is 0.965; the velocity 1.75 chord lengths both 
upstream and doTOstream is 0.290. 

APPEOXIMATE VELOCITY CORRECTIONS FOE AN AIEFOIL PLACED 
ALONG' CENTER LINE OF A PLANE-WALLED CHANNEL 

If an airfoil is placed midwaj' between the walls of a plane- 
walled channel, simple approximate velocity corrections may be 
derived under the conditions that the angle of attack is small and 
that the thickness, chord, and camber are small in comparison with 
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the dimensions of the channel. Tvo such corrections vill he 
exxDlained. Both corrections depend upon the successive reflection 
of the airfoil in the channel walls hy which a cascade of airfoils 
alternately upright and inverted is obtained. As is well knovni 
(see reference 9), the flow through such a cascade is equivalent 
to the flow about the airfoil in the plane-walled channel. In the 
first-order approximate theory, the image airfoils are replaced by 
doublets and by vortices; in the more elaborate treatment developed 
by Goldstein (reference 4), higher-order singularities are included. 
Inasmuch as the method of confomal mapping developed in the present 
paper is applied numerically to a symmetrical airfoil at the center 
of the channel, the approximte theories will be quantitatively 
discussed only for such airfoils. A more general treatment would 
follow along similar lines. 


First-Order Theory 


In the development of the first-order theory the vortex and 
the doublet are assumed to contribute independentlj'- to the velocity 
correction. The effect of the image vortices is to curve the stream, 
and to increase the effective angle of attack and lift on the airfoil 
in the charnel. The image doublets increase the velocity at the 
center of the channel and thus take into account the constricting 
effect of the channel walls. Glauert (reference 9, p. 49) obtained 
a fonnula for the ratio of the lift in the free stream to the lift 
in the channel. If it is assumed that the vortices merely change 
the angle of attack, the Kutta condition combined with Glauert 's 
formula yields the following result: 


sin a _ f £ f 

sin (x^ 24 [hj 


(42) 


where a is, as before, the angle of attack with respect to the 
direction of the flow at infinity and is the effective angle 

of attack due to the curved stream. 


The increase of velocity at the center of the channel induced 
by the image doublets is assumed to bo that due to the airfoil at 
its angle of zero lift. If this increase is denoted by u and, 
as before, T is the velocity at infinity in the channel, the 
following relation is true; 


u 

V 




( 43 ) 
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where for symmetrical airfoils 



as in reference 9 (p. 55). Here v^' is the velocity on the air- 
foil when the airfoil is in the channel at an angle of attack of 0° 
and y is the distance to the upper surface of the airfoii measured 
normally from the chord line. The integral in equation (44) is 
taken w5.th respect to the surface distance s along the upper sur- 
face of the airfoil from leading to trailing, edge. 

In the calculation of the strength of a doublet that is to 
replace an isolated airfoil, rather than ' should be used. 

However, inasmuch as the strength of the doublet must be increased 
when it is used to replace the same airfoil in cascade, the use of 
Vq', which is greater than Vj^, will change the value of x in the 
right direction. 


The velocity correction is defined as 



(45) 


where 



is the velocity on the airfoil in the channel expressed 


the airfoil is 
is the isolated 

airfoil velocity for the angle of attack a. Since the airfoil is 
small compared with the breadth of the channel, the flow about the 
airfoil in the channel is equivalent to the flow about an airfoil at 
an angle of attack in a free stream whose velocity at a great 

distance away is V + u. Therefore the following relation is true 


as a fraction of the ultimate upstream velocity when 
at an effective angle of attack and where 
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The fomula for the velocity correction shows the importance of the 
changed angle of attack, for one part of the correction is the 
difference in the isolated airfoil velocity distrihuions at angles 
of attack a and the other term of the correction is propor- 

tional to the isolated velocity distribution at the increased angle 
of attack. 


The correction obtained by the use of vortices and doublets 
is valid to the first order in and tc/h^. Wlaen the angle 

of attack is O'^, the pai’ameter j does not appear (reference l) . 


Goldstein's Second-Order Velocity Correction 

Goldstein (reference 4) first replaces the image airfoils by 
the doublet, the vortex, and the higher-oi’der singularities given by 
the potential function of the airfoil in a uniform free stream.. The 
nonuniform disturbance velocity produced by these singularities in 
the physical region, in particular at the location of the plvsical 
airfoil, is calculated, taking into account the change in direction 
of the stream. This first approximation nonuniform disturbance 
velocity (a) changes the velocity distribution on the airfoil from 
its isolated free-streara. value and (b) changes the value of the 
singularities that are to be imaged. Change (b) is evaluated and 
a second approximation nonuniform distribution of the airfoil in the 
final nonuniform, stream is calculated. 

In principle, Goldstein's method is capable of yielding to any 
degree of accuracy the effect of a plane-walled channel on the two- 
dimensional velocity distribution of an arbitrary airfoil, arbitrarily 
situated. The successive approximations become increasingly labo- 
rious, however, and only the second -approximation formulas are given 
in reference 4. 

The second -approximation formula for the constriction correc- 
tion for the symmetrical airfoil situated in the center of the 
channel at a small angle of attack is obtained as: 


\V/ 


/TTirP(e) - P(jt) + sin (0 + 


[ sin (9 + a) + 


tt]^) + sin aj 
sin a] 


(49) 


so that 



24 


NAG A EM No, E7A28 


Av 

V 



( 50 ) 


where U here represents the sum of the ultimate upstream velocity 
and the velocity at the center of the channel induced hy the singu- 
larities so that ^ 1 corresponds to u/v of the first-order 

theory; a-) is, as in the previous approximate theory, an effective 
angle of attack with respect to the direction of the stream; the 
function P(0) is a measure of the distortion of the stream caused 
by the singulax-ities. 


The Go3.dstein second-order image correction is accurate to the 


orders 


(tl 


! 7 


\hr\h)> v,4 ^ 


, 2.2 


ct' 

h^ 


and 1 - 


IIV 

Vh) 


VJhen the angle 


of attack is zero, the terms 
appear . 



do not 


Discussion of Numerical Eesults 


of Approximate Theories 


The first-oi-der and second-order corrections were calculated 
for the 12 -percent -thick sjunmetrical airfoil. The corrections for 
the airfoil at zero lift have been discussed in reference 1. The 
results for the angle of attack of 4*^ are plotted in figure 8. The 
constants used in the first-order correction are 



u/V 

°1 

3.S3J 

Q.0116 

4.458° 


Those for the second-order correction are 


^0 

Cl 

1 

1 

C2 

f 

Cs 

1 

C4 1 

- ... 1 

1 

i 

1 

1 

! V 

0.08722 

0.05534 

^0.02401 

0.00455 

0,00475 

4.200° 

0.0108 



NACA EM No. E7A28 


25 


The first-order theorj’’ yields good results for the upper 
surface of the airfoil in that the correction so derived shows the 
sanie over-all trend as the correction obtained by conformal mapping. 
The approximate correction appears to be a mean curve to which are 
added components due to the curvature of the airfoil. For the 
lower surface/ the approximate correction is not quite so good a 
mean line as it is for the upper surface. For both upper and lower 
surfaces, the contribution to the velocity correction due to the 
doublets and that due to the change in angle of attack are equally 
effective in forming the total correction. 

For the upper surface of the airfoil, the Goldstein second-order 
image correction follows the same trend as the first-order image 
correction, but the values are Diore nearly constant. The second- 
order correction for the lower surface follows more closely the trend 
of the mapping correction than the first-order correction. From this 
example,- the second-order correction appears to be more accurate than 
the first-order correction. 

The incremental velocities u/V and ^ - 1 of the first- and 

second -order corrections, respectively, are in good agreement but 
the values of the effective angles of attack a 3 _ differ markedly. 

This difference accounts for the difference in the nature of the 
correction curve of figure 8 near the leading edge of the airfoil. 


CALCULATION OF LIFT AND MOMEIff 

For the case of angle of attack of 4°, the lift coefficient c^ ' 
for the airfoil in the channel was calculated by integrating the 
pressure distribution about the airfoil. The calculation for C 2 ' 
was also carried out by means of the two approximate theories. 

The isolated airfoil lift coefficient C 2 was 0.478. The 
value of Cj ' obtained by the integration of the pressure distri- 
bution is 0.537; that value obtained from the second-order theory, 
0.522 by the formulas of reference 4; and that value obtained from 
reference 9 (p. 49), 0,532. All the values of c^ ' obtained indi- 
cate the expected increase in lift for the airfoil in the channel, 
and also show good agreement among themselves in that they do not 
vary more than 3 percent. The lift-coefficient correction, c^' - c^, 
varies, however, about 30 percent among the different theories. 

The lift coefficient C 2 ‘ was also calculated by integrating 
the pressure distribution on the walls of the channel. Theoretically, 
the integration should be carried out to infinity on either side of 
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the airfoil. The practical calculation is, of course, impossihle. 
The integration is therefore carried out only over a finite range 
to yield the lift coefficient Cj", and a correction factor used 
to take into account the effect of the rest of the channel. 

The correction factor t), which is equal to C2 "/cj*, has 
been derived in an approximate form, in the appendix of reference 10, 
The airfoil is replaced by a row of voi-tices, which are imaged in 
the walls of the channel. The- q factor for an individual vortex is 
calculated. The final t] factor is obtained by averaging q for 
each vortex with a loading derived from, thin airfoil theory as a 
weighting factor. 


In figure 9 the lift coefficient c^" is plotted as a function 
of the limits of integi'ation, wiiich were taken symmetrically about 
the origin. The value of c^”, obtained by integrating the pressure 
distribution 1.75 chord lengths upstream and downstream of the origin 
is 0,493. When this value is divided by the value of c^', derived 
by integrating the pressure distribution on tlio airfoil, a value t) 
of 0.918 is obtained. The value of q obtained by the method of 
reference 10 is 0.900. The value of c^', obtained from the approx- 
imate value of T], is 0,548. The correction factor obtained by the 
approximate method is satisfactory to the order of the approximate 
theories previously discussed. 

It is also possible to obtain the moment on the airfoil about 
any point by integrating the moment of the pressure (accurately 
calculated) on each element of area on the channel walls. A factor 
analogous to the rj factor can be so determined that the integration 
for the moment over a finite range may be extended to take into 
account the regions on the channel walls a great distance away. 


COriCLUSIOWS 

The analysis and numerical calculations of the present paper 
lead to the following conclusions: 

1. The method of conformal transformation by meeins of the 
Cartesian mapping function provides a satisfactory numerical solu- 
tion to the problem of obtaining the local velocity corrections for 
an arbiti’ary airfoil in a channel for the case of two-dimensional 
frictionless incompressible flow. 

2. If closeness to the velocity corrections obtained by con- 
foimal mapping is used as a criterion, the second-order Goldstein 
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correction is more accurate than the first-order image vortex and 
doublet correction for thin airfoils at small angles of attack in 
giving velocity corrections in the examples calculated. 

3. If it is necessary to obtain a higher-order correction than 
the second^ the method of the Cartesian mapping function is probably 
more convenient to use than the Goldstein type correction. 

4. The channel lift coefficients obtained by the two approxi- 
mate theories are in good aigreement with the lift obtained from the 
mapping velocity distrib\ition; the lift corrections obtained by the 
two approximate theories are not in good agreement with the correc- 
tion obtained by mapping results. 

5. The existing method of finding the lift coefficients from 
the velocity distribution on the channel walls has been satisfac- 
torily checked. 


Aircraft Engine Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, December 4, 1946, 
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APPEiroiX A 

DEEIVATION OF THE EELATIONS BETWEEN THE PEAL 
MD IMAGIWAEY PAETS OF TEE aiF 

InasEii’ch as the C1’'IF z - ^ is regular within the annulus and 

also on the hounding circles in the p-plane^ it xnay he expanded in 
a Laurent series, which is valid in the axinulus and on the circles 
hounding the annulus. Thus 



J? 

z - C = Ax(p,©) + iA 7 (p,-p) = (Al) 

From eciuation (Al) the following expressions are obtained: 

00 

AX]_(c?) = ao (a^ + a_^^) cos ncp- (h^, - h_r^) sin ncp 


00 00 


Ayi(9) =tQ -;-2l_.(^n -s-n) sin n.p+ ^'(hn + h_n) cos ncp 

1 ' 


1 

OO 


? (A2) 


Ax2(9) = Sq +y ~ !> " ^-n*!”") sinnp 


1 

OO 


1 

OO 


Ay 2 (c?) =to a.^q"^) sin n^ +N + h.^q"”) cos ncp 

1 ' ”T“' 




The values of a^^ and h^^ can he found hy raeans of Fourier's rule 
in tenas of the CMF. 

Vfiien aQ and hQ are evaluated, the conditions of consistency 
that are necessary conditions for the regularity of the CMF in the 
annulus appear as 

f 2it f‘ 2n 

Ax^ (cp) dy = I AX2('-P) d:p 


2naQ = 


JO 


2rthQ = 


J 


■27T 

Ayi (cp) = 

o'" ^ 


2it 


AygCq?) <19 


(A3) 
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Four relations are desired : and AXg expressed in terras 

of Ay-]^ and Ajg and, conversely, Ay]_ and expressed in 

temis of Ax^ and Axg. The derivation of the expression for Ayi^ 
in terms of Ax^ and Ax£ will now he carried out. The other 
relations vrill follow analogously. 


Through the use of the first and third equations of (A2) and 
through the use of Fourier's rule, the coefficients a^^ and hj^ are 
evaluated. As a result of the calculation, the following equations 
are obta j ned : 


-D^q'^ + Dg 
qn _ g-n 


a 


-n 


- Dg 1 

qn . q^-n j 


n 


- % 

qn - q-n 


h 


-n 


qn - q“n 


(A4) 


where 




D 


2 


r 2rt r 2rt 


1 ' 
.'0 

AXq(-p) 

cos ncp dcp 

-Jo 

f2n 



r2rt 

1 ( 
Jo 

AxgCcp) 

cos n9 d<p 

Kg = - 

^ ^Jo 


AxT^(cp) sin ncpdcp 


> (A5) 

' 

AX 2 (cp) 3in ncpdcp 


The values of the coefficients and are substituted in the 

infinite series expression fox’ given by equations (A2). The 

values for ^2 given hr equations (AS) are also 

used. Thus 


Ay,(cp)=hQ+— y — j Ax-^( <p') (sin ncp' cos ncp - cos n cp' sin ncp)dcp' 

qn- q-njo 


K 


\ 


, - 1' 


•n 


2jc 


Ax 2 ('P') (cos ncp' sin n<p - sin ncp' cosncp)dcp' 


_/ 


1 
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00 

1 \ n ~n r 

Ayi(cp) = '^0 + ^ y ^n"~^- ' n jQ Ax^^Ccp') sin n(q>‘-cp) dcp« 


CD 

r V 


1 \ ' ? 

- i ) — / Ax^Ccp') Sin n(cp'-cp) dcp' 


(A6) 


Now let f (cp) te a function that can be developed in a Fourier 
series for 0 ^ cp< 2jt. Then 



.V 

■/, 

m=l 


P 


f (cp> ) sin n(cp‘ -cp) 


dcp' = 0 


(A7) 
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Inasmuch as the series of equation (A8) are uniformly convergent as 
are the series of equations (5), the summation and integration may 
he interchanged in equation (A8) to yield equation (5) . 
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APPENDIX B 


TEE NUMERICAL EVALUATION OF THE CARTESIAN MAPPING FLUJCTION 


The determination of the functions AX]_ and AX 2 from the 
given functions and Ay 2 was "based in this paper on nuiaerical 

integration of the first two equations (5). The equations for AX]_ 
and AXg^ when the constant aQ has "been set equal to zero, are 


Ax^(cp) = i 


r 2n 


Ayi( <?') 


1 cot .Y Y- 

2 2 , i_q2n 

1 


d(f ’ 


2jt '’s 

i I Ayo(9')N — sin n(9’-9) d9' 

-Jo ^ / . l-pg-i 


r 


AXgCcp) = - i 


2jt 


AygCcp’) 


“1 


isothel# YYi n(9'-9)|a»' 


Z - 1-q 

1 


2n 


+ - Ay-,(cp*)\ — ^ sin n(cp'-q5) dcp' 

i 


If the range of cp is divided into 2n equal intervals whose length 
is 5, if the values of Ay are given at the end points of the 
intervals, and if Ax is desired at the same points, approximate 
integration will yield expressions of the following form: 


Ax^(cp) 


Ax^Ccp) 


2n-l 


2n-l 




[cj^ + djj.) Ayj_(cp + k5) + 

■■ 

®k 

k5) 

it 

O 


k=0 



2n-l 


2n-l 



'V ^ 

. 

(cjij. +.di^) AygCcpd- k5) - 
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/ 

^Ylicp + 

k5) 

k=0 


k=0 




> (B2) 
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The values of C]j have been calculated in reference 1 by means of 
Simpson's rule and other simplifications for use with the CMF of 
simply connected regions. The values of di, and e> may be 

similarly obtained. The value of cj^ as calculated in reference 1 
is 


cq 


= 0 


S 6 S + sin 5 

° 2 2n sin 6 


2n-l 6jt 


-S b _ 5 -I- sin S 


2n sin 5 


(B3) 




f cotf 
37t 2 


(k odd) 


Cj, = 


26 , k5 

— cot 
3rt 2 


(k even) 


I 

In the present paper, because the number of intervals was an integral 
multiple of 6, Weddle's rule was used for the evaluation of dj^- 
and ejj . 


The values of cjj. are given in table 4 foi- the cases of 
2n = 24 and 2n = 48. The values of and ejj. contain the 

paraiaeter q,. Hence, these coefficients must be evaluated anew for 
each approximation. 
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TABLE - ORDII'IATES OF 12-PEECEMT THICK AIRFOIL 
[ From table i of reference l] 


Station 

(percent 

chord 

from 

nose) 

Ordinate 

Station 
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chord 
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nose) 

Ordinate 

. 

0 

0 

50 

5.880 

1,25 

1.425 

55 
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2.5 

1.900 
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5.025 

5 
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65 

4.415 

10 

3.540 

70 

3.750 

15 

4.250 

75 

3.060 

20 

4.820 

80 

2.350 

25 

5.295 

85 

1.685 

30 

5.655 

90 

1.060 

35 

5.900 

95 

.510 

40 

6.000 

97.5 

,260 

45 

6.010 

100 

i 
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TABLE 2,- VELOCITY DISTRIBUTION AND CARTESIAN- MAPPING FUNCTION FOR AIRFOIL AT ANGLE OP ATTACK OP 0 
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TABLE 3.- VELOCITY DISTRIBUTION AND CARTESIAN MAPPING 
FUNCTION FOR AIRFOIL AT ANGLE OP ATTACK OP 4® 
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TABLE 4. - COEEFICIEIirS FOR CALCULATION OF CARTESIAN 

MAPPING FUNCTION FOR SINGLE CONTOUR 
(a) 24-point scheme 
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^k 
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°k 
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12 

0 
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.42564 

13 

-.00366 
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.20734 

14 

-.01489 

kJ 
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(h) 48-point 

scheme 
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Figure 1,- Transformation of flat plate and 
channel Into two conoentrlo circles. 



Channel wall In z-plane 



gure 2,- Transformation of airfoil and arbitrary channel Into flat plate 
and plane-walled channel by Cartesian mapping function. 




Channel 







Velocity dlstrlbutlonf ^ 


Fig. 4 
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Figure 4,- Comparison of velocity distributions In channel 
obtained by two methods. 
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Fig. 5 
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Figure 5*- Velocity distribution on airfoil In channel at angle of attack of 4 ® and velocity distribution on channel walls*. 



Velocity distribution* Vj^/V. 
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Figure 6,- Comparison of velocities on airfoil In free 
stream and on airfoil In channel for a > 4°. 




Fig. 7 
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Figure 7,- Cartesian mapping function for airfoil In a 
channel at 'c » 4<>. 



Local velocity correction, Av/v 


Fig. 8 
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(b) Lower surface* 

Figure 8,- Velocity corrections for 12-percent-thlck airfoil. 
0=4°, c/h. = 0,5. , 
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Fig. 
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